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Abstract 

Collapse and fragmentation of uniform filamentary clouds under isotropic far- 
ultraviolet external radiation are investigated. Especially, impact of photodissocia- 
rh \ tion of hydrogen molecules during collapse is investigated. Dynamical and thermal 

evolution of collapsing filamentary clouds are calculated by solving virial equation 
Ph. and energy equation with taking into accounts non-equilibrium chemical reactions. 

It is found that thermal evolution is hardly affected by the external radiation if the 
initial density is high (no > 10 2 cm -3 ). On the other hand, if line mass of the filamen- 
tary cloud is moderate and initial density is low (no < 10 2 cm -3 ), thermal evolution 
of the filamentary cloud tends to be adiabatic owing to the effect of the external 
,_j- \ dissociation radiation. In this case, collapse of the filamentary cloud is suppressed 

and the filamentary cloud fragments into very massive clouds (~ 10 4_5 M Q ) in the 
^sO ■ early stage of collapse. Analytic criterion for the filamentary clouds to fragment into 

such massive clouds is discussed. We also investigate collapse and fragmentation of 
CN ' the filamentary clouds with an improved model. This model can partly capture the 

effect of run-away collapse. Also in this model, the filamentary clouds with low initial 
density (no < 10 2 cm -3 ) fragment into massive clouds (~ 1O 4 M ) owing to the effect 
;_! \ of the external radiation. 
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1. Introduction 

It is widely accepted that in the dark age the first collapsing density perturbation col- 
lapses and cools owing to hydrogen molecules (H2) to form so-called population III (popIII) 
stars (Bromm et al. 1999, 2002; Abel et al. 2000, 2002; Yoshida et al. 2008). PopIII is expected 
to form in halos typically with ~ 10 6 M Q (Tegmark et al. 1997). If popIIIs are massive stars, 



they are expected to cause feedbacks via supernovae and radiation. The former may spread 
metals made by nuclear fusion and sweeps neighboring gas by shock (e.g., Heger et al. 2003). 
This paper focuses on the latter. There are two types of the radiative feedbacks, which are 
ionization and dissociation (Whalen et al. 2004; Susa et al. 2009; Hasegawa et al. 2009a, 
2009b). Ionization of hydrogen atoms (H) provides strong heating which causes evaporation of 
clouds. However, since inter stellar matter mainly consists of H with large opacity, ionization 
photon tends to be prevented from spreading, and photoionization would occur mainly within 
halos. On the other hand, photodissociation of H2 would occur even out of halos (Kitayama 
et al. 2004). Thus, it is expected that there are some regions which are not photoionized but 
photodissociated. We investigate gravitational collapse and fragmentation of primordial clouds 
in such a region. 

Since H 2 is main coolant in the early universe, if H 2 is photodissociated, a collapsing 
primordial cloud (> 1O 8 M ) heats adiabatically up to high temperature (~ 10 4 K) where atomic 
cooling becomes effective. In such a case, it is suggested that high mass objects (~ 1O 4_5 M ) 
may form because Jeans mass [Mj oc T 3 / 2 p -1 / 2 ) is large (e.g., Bromm & Loeb 2003). This 
possibility may connect the direct formation of supermassive black hole. However, if collapsing 
cloud fragments into many small clumps, collapsed object may form star cluster (Omukai & 
Yoshii 2003). In order to understand actual final outcome and their initial mass function 
clearly, further detailed investigation of thermal and dynamical evolution of collapsing clouds 
is required. 

Omukai (2001) investigated the evolution of spherical clouds under the dissociation radi- 
ation. The author calculated the evolution of the central region of a collapsing sphere, assuming 
free-fall collapse. Started from no = 8.9 x 10~ 2 cm _1 , the clouds collapse adiabatically in the 
early stage, n < 10 2 cm -3 . After this stage, thermal evolution of clouds is divided into two 
types of tracks. When the external radiation is very strong (e.g, with the intensity larger than 
10~ 18 erg cm _2 s _1 Hz _1 sr _1 at 13.6eV for thermal radiation of 10 4 K), H 2 is photodissociated 
enough to suppress H 2 cooling. In this case, main coolant is provided by hydrogen atom. On 
the other hand, if the intensity of the external radiation is moderate, sufficient amount of H 2 
forms and clouds cool mainly via H 2 cooling. Susa (2007) investigated more realistic evolution 
of a spherical cloud under the UV radiation from a single light source by three-dimensional 
calculations. The author investigated whether or not clouds collapse for parameters such as 
distance from the light source and the density when the light source turns on. 

As for the formation of the spherical clouds investigated above, a filamentary cloud is a 
possible origin. Filamentary clouds are commonly expected during the way to form the stars. 
When sheet-like cloud forms, the sheet-like cloud tends to fragment into the filamentary clouds 
(Miyama, Narita, & Hayashi 1987a, b). In cosmological simulation, the filamentary structure 
forms from density perturbation which has > 10 6 M Q (Abel et al. 1998; Bromm et al. 1999; 
Greif et al. 2008). These filamentary clouds have a possibility to produce spherical clouds by 



fragmentation (Nagasawa 1987; Inutsuka & Miyama 1997). Thus, in order to understand the 
origin and property of a collapsing spherical cloud, it is useful to investigate the evolution of fil- 
amentary clouds. Previous works about fragmentation of primordial filamentary clouds include 
one- zone models (Uehara et al. 1996; Flower 2002; Omukai & Yoshii 2003) and one- dimensional 
models (Nakamura & Umemura 1999, 2001, 2002; Uehara & Inutsuka 2000). Among these stud- 
ies, Uehara et al. (1996), Omukai & Yoshii (2003), and Nakamura & Umemura (1999, 2001) 
considered only H_2 as a coolant, and the others considered deuterated hydrogen molecules (HD) 
as well as H_2. If the initial fraction of H2 is lower than 10 -3 , H2 becomes main coolant and 
the fragment mass is 1 — 5OOM . On the other hand, if the initial fraction of H 2 is higher 
than 3 x 10 -3 , a filamentary cloud with low initial density (uq < 10 4 cm -3 ) cools mainly via HD 
cooling since H2 promotes the formation of HD with larger cooling rate. In this case, the frag- 
ment mass is 1 — 140M Q . When HD becomes main coolant, the filamentary clouds cool down 
to lower temperature (e.g., ~ 40K) and fragment into less massive fragments 1 than when H2 
is main coolant. All of these previous papers except Omukai & Yoshii (2003) did not consider 
the feedback effect from the external radiation. 

Once a massive star forms in a cloud ~ 1O 6 M , the whole cloud is photodissociated by 
UV radiation from the star (Omukai & Nishi 1999). In order to form subsequent stars, larger 
cloud e.g., > 10 8 M Q is required. Omukai & Yoshii (2003) investigated fragmentation of the 
filamentary cloud under the UV radiation using a one-zone model. Based on the results by 
Nakamura & Umemura (2001), they assumed that the filamentary cloud fragments when its 
density becomes 100 times higher than that at the loitering point at which temperature is a local 
minimum in p — T plane owing to H 2 cooling. The authors concluded that the fragment mass 
is smaller under the stronger UV radiation. However, the condition for fragmentation given by 
Nakamura & Umemura (2001) is considered for the cases without the external radiation. Thus, 
it is not clear that this condition for fragmentation is applicable for the cases with the external 
radiation. When the filamentary cloud suffers photodissociation, it will collapse adiabatically 
and it will fragment before it reaches loitering point. In such a case, the fragment mass becomes 
larger than in the cases without the external radiation. 

In this paper, we investigate the effect of photodissociation radiation on the thermal 
evolution of the collapsing filamentary cloud and on the mass of fragment. We investigate 
whether or not the external radiation photodissociates enough H2 to prevent H2 cooling and 
whether or not the external radiation makes fragment mass larger. Instead of assuming free-fall 
collapse, we calculate the dynamical evolution by solving virial equation taking into account 
the effect of pressure gradient force as a result of insufficient cooling due to photodissociation. 
As for the condition for fragmentation, we assume that the filamentary clouds fragment when 
the timescale of fragmentation becomes shorter than the timescale of density evolution (see 



1 Since low temperature helps collapse, the timescale of density evolution is short and fragmentation does not 
occur until density becomes very high (~ 10 12 cm -3 ) (see §2.3). 
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§2.3). 

We measure the effect of the external radiation with respect to fragment mass. We 
investigate the dependence of fragment mass on parameters such as initial density, line mass, 
and intensity of the external radiation. We also discuss analytically what physical process 
determines fragment mass. We consider the simplest case in which the dissociation radiation 
is isotropic. As for the intensity of the external radiation originating from forming popIIIs, we 
refer the result by Dijkstra et al. (2008) which investigated the mean intensity of the external 
dissociation radiation in the universe at z ~ 10. 

In §2, we describe the model which is used in this paper. We present numerical and 
analytic results in §3 and §4, respectively. We present an improved model and numerical 
results in §5. Finally, §6 is devoted to summary and discussion, including discussion about 
evolution of cloud after fragmentation. 

2. Model 

2.1. Basic equations 

For simplicity, we assume that the filamentary clouds are uniform. We solve the virial 
equation for the dynamical evolution in the cylindrical radial direction (Uehara et al. 1996). 
We do not consider dark matter for simplicity. In the case with high initial density, baryon 
density is expected to dominate dark matter density and it will give a good approximation. In 
the case with the low initial density, we will underestimate the effect of dark matter, e.g., large 
infall velocity owing to dark matter gravity. The virial equation for the filamentary cloud of 
unit length with volume V is 

where G is gravitational constant, / is the line mass (mass per unit length) of the filamentary 
cloud, 

I=j pr 2 dV (2) 

is the inertial moment per unit length with density p and radius r, 

* = /Jp«W (3) 

is the kinetic energy per unit length with the velocity v, and 

n = f PalV (4) 

J \ 

is the integrated pressure per unit length with local pressure P. By substituting equations 
(2)-(4) into equation (1), we have 
d 2 R 2G 



alt 2 R 



{l-lc(T)h (5) 



where R is the radius of the filamentary cloud, T is temperature, and l c {T) is the critical line 
mass for the hydro-static isothermal filamentary cloud defined as 

« r > - IS (6) 

(Ostriker 1964) with Boltzmann constant ks, mean molecular weight p, and mass of a hydrogen 

atom 77tH- 

For the thermal evolution, we solve energy equation 

du _ _ d 1 A radi A chem 

dt dt p p p 

where u is thermal energy per unit mass 

1 k B T 

u = 8 

7 ad - 1 pm H 

with adiabatic exponent 7 a d- The first term of the right hand side in equation (7) denotes 

adiabatic heating. The radiative cooling rate A ra di per unit volume includes lines of H, lines of 

H 2 , lines of HD, and continuum (see tablel in detail; see also Omukai 2001). For the radiative 

cooling, the effect of radiative transfer is included according to Susa et al. (1996). The symbol 

Achem represents heating/cooling rate associated with chemical reactions. Equation of state for 

ideal gas 

P=PML (9) 

pm u 

is assumed. 

We consider non-equilibrium chemical reactions by solving equations 

-jT = Y, k ijkfjfkn + Y, hjfj, (10) 

al j, k j 

where n is number density of all nuclei, kijk and fcy are reaction rates of formation and destruc- 
tion of species i, and fi is the fraction of species i. We consider the following fourteen species : 
H, H+ H", H 2 , Et, He, He+, He ++ , D, D+, D~, HD, HD+, and e~. We consider 35 reactions 
concerned with H and He taken from Omukai (2001) and 18 reactions concerned with H and 
D taken from Nakamura & Umemura (2002). We also consider photodestruction of species D, 
D - , and HD+ (we refer to Galli & Palla 1998 for D and HD+ and Frolov 2004 for D~). Above 
equations are solved numerically with implicit integrator. 

2.2. External radiation 

We assume the external radiation to be isotropic. Dijkstra et al. (2008) investigated the 
mean intensity of the dissociation radiation at z ~ 10 from the surrounding star-forming halos 
and estimated the probability distribution of the mean intensity. We adopt mean intensities 
whose probabilities are ~ 0.4 2 and ~ 0.06 in Dijkstra et al. (2008) (see §2.4). Moreover 



0.4 is the highest probability. 



we assume that the external radiation is thermal radiation from 12OM stars and surface 
temperature (T sur = 95719K) of stars is determined according to Schaerer (2002). We also 
assume that the ionization photon does not reach the filamentary clouds from light sources. 
When we calculate photodissociation reaction 

H 2 + 7 ^H 2 *^2H (11) 

(solomon process) 3 , we consider the extinction of photon by dissociation of H2 and absorption 
by continuum processes (see Table. 1). The photodissociation rate is proportional to the mean 
intensity of the dissociation radiation. During penetrating the filamentary cloud from the 
surface to the center, the intensity of the dissociation radiation decreases owing to dissociation 
of H2 and absorption by continuum processes. In this paper, the effects of radiative transfer 
of dissociation photon is approximated by the product of shield factor as J v = f a hfconJv,o, 
where f s ^ is self-shielding factor associated with photodissociation of H 2 , f con is decreasing 
rate associated with absorption by continuum processes, and J Vt o is the mean intensity of the 
dissociation radiation at the surface of the filamentary cloud. 

First, we consider the photon decreasing rate, / con , associated with the absorption of 
dissociation photon by continuum processes. We focus on the dissociation photons with 12.4eV. 
Radiative transfer equation along the s-direction is given by 

^ = -W„ + j„, (12) 

where I v is the intensity of the radiation of frequency u, k u is the total opacity associated with 
reactions in Table. 1 (see Appendix 1), and j u is emissivity. For simplicity we assume that the 
scattered photons are absorbed immediately. We also assume j u = since j u mainly consists of 
lines and continuum with lower energy than 12.4eV. We consider the length of the column in 
various directions. Using the length R/sin8 from surface to the center of the filamentary cloud 
with angle 9 from the axis of the filamentary cloud, the intensity of the external radiation at 
the center is given by 

I u (0) = J ufi exJ-^-). (13) 

V sin fly 

Hence, f con J vfi is given by 

fconJvfl = -jr- d(j) d9 sin6»exp( — ?—). (14) 

Air Jo Jo V sin 6'/ 

The integral of the right hand side in equation (14) is calculated with the fitted functions given 
in Appendix 2. In typical results shown in §3, the value of / con is found to be larger than 0.97 
during collapse for n < 10 2 cm~ 3 , where the photodissociation of H 2 is effective. The absorption 
in the low density cloud has only a minor effect. 

Second, we consider self-shielding factor, / s h, associated with the dissociation. The 
self-shielding factor / sh is approximated by 

3 H2* is excited H2. 7 is photon with 12.4eV. 



/ sh = min 



A^h 2 \" 3 / 4 



(15) 



10 14 cm- 2 , 

(Draine & Bertoldi 1996), where Nn 2 is the column density of H 2 . To estimate the effective 
column density, we estimate average in angle of the length between the surface and the center. 
Using the length R/sin8 from the surface to the center of the filamentary cloud in the direction 
with angle 9 from the axis of the filamentary cloud, effective column density of the filamentary 
cloud is estimated as 

\ i-2-tr r-K R TV 

7V H = t~ / d(j) d8 sm6^—n H2 = -n K2 R, (16) 

An Jo Jo sin u 2 

where nu 2 is number density of H_2. Since R oc n" 1 / 2 , effective column density increases during 

collapse as Nn 2 oc rtu 2 R oc n 1 / 2 . Finally, the photodissociation reaction rate of H 2 is given by 

k 2step = 1.4 X lOVcon/sh^O S" 1 . (17) 

2.3. Fragmentation of filamentary clouds 

During collapse of the filamentary cloud, two important timescales exist. One is the 
timescale of density evolution defined as t& yn — Pel ' Pc where p c is the density at the center. The 
other is the timescale of fragmentation defined as £f rag = 2.l/y/2irGp c (Nagasawa 1987; Uehara 
et al. 1996). The latter is the timescale in which the fastest growing mode of perturbation 
grows to non-linear. According to Inutsuka & Miyama (1992), if acceleration in the radial 
direction is high, perturbation with low amplitude does not grow enough during collapse. When 
acceleration in the radial direction becomes low owing to strong pressure, tdyn becomes large 
and the growth of perturbation becomes remarkable. The wave length of the fastest growing 
mode of perturbation is nearly the diameter of the filamentary clouds (Nagasawa 1987). This 
wave length becomes shorter during collapse. If the fastest growing mode has time to grow to 
non-linear before the diameter of the filamentary cloud changes largely, the filamentary cloud 
is expected to fragment. Thus, we assume that the filamentary clouds start to fragment at the 
moment when £f rag < tdyn is satisfied (Uehara et al. 1996; Inutsuka & Miyama 1997). 

Using the wave length of the fastest growing mode Af rag ~ 2ixR at fragmentation, the 
fragment mass is estimated as 

M frag = A frag i ~ 2ttRI (18) 

(Narita, Miyama, & Hayashi 1987a, b; Larson 1985; Uehara et al. 1996). According to equation 
(18), the fragment mass is proportional to the radius of the filamentary clouds. If fragmentation 
occurs after the filamentary cloud collapses to a small radius with high density, the fragment 
mass is small. 

2.4. Parameters and initial conditions 

In this paper, we treat three physical quantities as parameters, which are initial number 
density n , normalized intensity of the external radiation, 
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T _ Jhv=13.6eV,0 , lrA 

J<21 ~ in 91 9 VLJ 1 i"' l iy J 

and the line mass parameter, 

_ 7rGp fim n 2 / , , 

; = -^^^^ " o^y (20) 

with initial density p ; initial temperature T , and initial radius Rq 4 . The reason why we 
choose these three quantities is as follows: in the view point of dynamical evolution, the line 
mass parameter / is important. In the view point of thermal evolution, initial density uq is 
important. The symbol J Vj q and n are necessary to study the effect of dissociation photon. 

We consider cases with log 10 ra = 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5, 5.5, and 6 for n and 
/ = 1.25, 1.5, 1.75, 2, 2.25, 2.5, 2.75, and 3 for /. For J 21 we consider J 21 = 1, 6.5, and 10. The 
case with J 2 i = 1 demonstrates the weak external radiation. According to Dijkstra et al.(2008), 
J 2 i = 6.5 is the average intensity at z ~ 10, and J 2 i = 10 represents strong radiation case whose 
probability is 0.06 (see §2.2). 

We assume that radial infall velocity at the surface of the filamentary cloud equals to 
the sound speed. If the filamentary clouds form from the sheet-like cloud, gravitational force 
dominates pressure gradient force in the filamentary cloud. Thus, the infall velocity when the 
filamentary cloud forms is expected to be in the same order as the sound speed, i.e., v = ac s with 
a numerical coefficient a ~ o(l) which depends on the details of fragmentation. As a typical 
value, we set a = 1 according to Nakamura & Umemura (2002). We calculated several cases 
with various values of a and found that the evolution hardly changes for a < 5. We mention 
the case with a = 5 in §3.4. 

In this paper, the filamentary cloud is assumed to form from a cloud which experiences 
H 2 cooling without UV radiation. We also assume that the external radiation turns on when 
the filamentary clouds form. The initial values of temperature and /h 2 are set to be 300K and 
/h 2 = 10~ 4 . In addition to /, n , and J 2 i, for the thermal evolution, T and fraction of H 2 , /h 25 
are also important. We adopt one typical value for them. This value of fn 2 is typically seen 
in cosmological simulations (e.g., Abel et al. 1998), /h 2 ~ 10~ 4 — 10~ 3 . We discuss initial H 2 
fraction with the effect of the dissociation radiation in §3.5. Fraction of He, /fie, is set to be 
0.0825 which corresponds to the mass fraction Y p = 0.244 (Izotov & Thuan 1998). Initial fraction 
of electron, f e , is set to be 10~ 4 according to Uehara et al. (1996). We adopt this value of f e 
in order for electron not to change /h 2 largely in the early stage of collapse 5 . Initial fraction 
of proton is determined from the charge conservation. We assume that [D]/[H] = 4 x 10 -5 , 
which is consistent with observations of the deuterium Lya feature in the absorption spectra 
of high-redshift quasars (e.g., O'Meara et al. 2001). Fraction of the other species is set to be 



4 The typical value of / is 2. This value is realized when the sheet-like gas fragments by the fastest growth 
rate (Miyama et al. 1987a). 

5 Electron helps H2 formation via H~ channel. 
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zero at the initial state. 
3. Results 

3.1. Cases without the external radiation 

3.1.1. Low density filamentary clouds with large line mass 

Figure 1 shows the result of the case with low initial density and large line mass, 
(/, n , J21) = (3, 10cm -3 , 0). From the early stage of collapse, adiabatic heating rate is a little 
higher than H 2 cooling rate and temperature gradually increases. At n ~ 10 8 cm~ 3 , l c overcomes 
/, and the filamentary cloud begins to be decelerated. Above n ~ 10 8 cm -3 , the three-body re- 
action becomes efficient and H 2 fraction increases to ~ 0.4. At n ~ 10 10 cm -3 , although H 2 
cooling is still effective, chemical heating associated with the three-body reaction of H 2 for- 
mation also becomes effective. Thus, temperature continues to increase. When the density 
reaches n ~ 10 11 cm~ 3 , although the three-body reaction is inefficient, temperature stops to 
increase owing to sufficient cooling with a large fraction of H 2 . The filamentary cloud becomes 
optically thick to H 2 line emissions at n ~ 10 12 cm -3 . Around this density, temperature in- 
creases again and it eventually exceeds 2000K. Such a high temperature state causes collisional 
dissociation of H 2 . Since chemical cooling associated with this dissociation can not dominate 
adiabatic heating, temperature is kept high enough to decelerate collapse. As a result, the 
filamentary cloud fragments when density reaches n ~ 10 15 cm -3 . Several authors pointed out 
that H 2 collision-induced emission becomes effective at n ~ 10 15 cm -3 (Omukai & Nishi 1998; 
Ripamonti & Abel 2004; Yoshida et al. 2006). However, in the case in figure 1, since tem- 
perature is high (~ 3000K), 80% of H 2 is dissociated and cooling rate of H 2 collision-induced 
emission is smaller than adiabatic heating rate by two orders of magnitude. Since the density 
of the filamentary cloud at fragmentation is very high (n ~ 10 15 cm -3 ), the mass of fragment is 
small (~O.1M ). 

In summary, the evolution of the low density models with large line mass is affected 
largely by radiative cooling and chemical heating/cooling associated with H 2 . In this sense, 
our result is qualitatively same as the previous results by Uehara et al. (1996) and Nakamura 
& Umemura (1999, 2001, 2002). Note that the above result with sub-solar mass of fragment 
originates from the one- zone model with a uniform filamentary cloud. In §5, we show the result 
with an improved model with the effect of run-away collapse. 

3.1.2. High density filamentary clouds with small line mass 

Figure 2 shows the result of the case with high initial density and small line mass, 
(f,n ,J 2 i) = (1.25, 10 6 cm~ 3 ,0). In this case, adiabatic heating dominates cooling a little after 
the early stage of the collapse, n <3 x 10 6 cm -3 . Collapse is accelerated only in the early stage 
of collapse (n < 2 x 10 6 cm~ 3 ) and not after that. Since acceleration is limited in the short 
density range, collapse of the filamentary cloud is limited at lower density {n ~ 10 8 cm^ 3 ) and 
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the fragment mass is larger (~ 50M Q ) than the case in figure 1 (§3.1.1). Different from the case 
in figure 1, H2 cooling never dominates. 

Figure 3 shows the fragment mass for various no and /. All lines are similar to each other 
and can be approximated as Mf rag ~ 230n ~ ' 03 / -5,1 with an error at most factor 4 at / = 3. This 
approximate function agrees with numerical results at low / (< 2). The fragment mass is deter- 
mined mainly by /. This tendency agrees with the result of Uehara et al. (1996). Nakamura & 
Umemura (2002) suggested that the fragment mass depends mainly on no- However, our results 
do not agree with that of Nakamura & Umemura (2002). This difference comes from simplicity 
that the filamentary clouds is assumed to be uniform. In the uniform model, virial temperature 
is determined by the whole line mass (/(= /// c )) 6 - The evolution of the non-uniform filamen- 
tary cloud (e.g., the one- dimensional model) includes run-away characteristics of the flow. The 
improved model with the effect of run-away collapse will be introduced in §5. 

3.2. Cases with the external radiation 

3.2.1. Low density filamentary clouds with large line mass 

Figure 4 shows the result of the case with low initial density, large line mass, and strong 
external radiation, (/,n , J 21) = (3, 10cm~ 3 , 10). The case in figure 4 corresponds to the case in 
figure 1 (§3.1.1) with the external radiation. In figure 4, it is seen that /h 2 decreases owing to 
photodissociation in the early stage of collapse, and that adiabatic heating dominates from the 
early stage of collapse. Cylindrical collapse is decelerated at n ~ 10 2 cm~ 3 since temperature 
increases. However, this deceleration is temporary and the filamentary cloud does not fragment 
at this point. Instead, it continues to collapse and shields itself from the dissociation photon. 
Then /h 2 begins to increase at n ~ 10 2 cm~ 3 . After that, H2 cooling becomes efficient and 
the evolution becomes similar to that in figure 1 (§3.1.1). As a result, the filamentary cloud 
collapses until it becomes optically thick to H 2 lines, and it fragments into the low mass clumps 
about O.14M . This mass of fragments is expected to be underestimated owing to the uniform 
filament model with homologous collapse as in the case in figure 1 (§3.1.1). 

3.2.2. Low density filamentary clouds with small line mass 

Figure 5 shows the result of the case with low initial density, small line mass, and 
strong external radiation, {f,no,J 2 i) = (1.25, 10cm -3 , 10). In this case, the external radiation 
photodissociates H2 in the early phase since dissociation photon penetrates the filamentary 
cloud with low column density. The early photodissociation suppresses H 2 cooling. As a 
result, temperature increases adiabatically until fragmentation. Since collapse is terminated and 
fragmentation occurs at low density (~ 34cm~ 3 ), the fragment mass is very large (~ 10 5 M Q ). 



(1 



Since collapse of the uniform filamentary cloud is homologous, virial temperature is determined by whole 
line mass, that is /. On the other hand, since collapse of the filamentary cloud is run-away collapse in 
one-dimensional model, virial temperature is determined by mass of the central region. The mass of the 
central region mainly depends on no- 
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The difference between cases with and without the external radiation is whether or not H2 is 
dissociated by the external radiation in the early stage. If H2 is photodissociated sufficiently 
enough to suppress H2 cooling, the filamentary clouds with low / (< 2.5) evolve adiabatically 
and fragment into very massive clumps. 

The difference between figure 5 and figure 4 (§3.2.1) is the value of line mass. Since 
the initial density is low, in both cases H2 is dissociated and the filamentary clouds evolve 
adiabatically in their initial stage. Since the line mass for the case in figure 5 is smaller, the 
slight increase of temperature is sufficient to suppress collapse. On the other hand, since the 
line mass for the case in figure 4 is larger, the slight increase of temperature is not sufficient to 
suppress collapse. Hence in the case in figure 4 the filamentary cloud does not fragment in the 
early stage of collapse and eventually H2 forms enough to cool the filamentary cloud. Difference 
of these two results originates from the line mass of the filamentary clouds. A critical line mass 
to shield themselves from the dissociation photon is discussed analytically in §4.4. 
3.2.3. High density filamentary clouds with small line mass 

Figure 6 shows the result of the case with high initial density, small line mass, and strong 
external radiation, (f,no, J21) = (1.25,10 6 cm _3 ,10). The case in figure 6 corresponds to the case 
in figure 2 (§3.1.2) with the external radiation. The evolution of the filamentary cloud in figure 
6 is similar to that in figure 2 (§3.1.2). This is because the initial density is high enough to shield 
the filamentary cloud from dissociation photon. Adiabatic heating dominates H 2 cooling from 
the early stage of collapse, and temperature increases gradually. In this case, the filamentary 
cloud fragments into the slightly more massive fragments than in the case without the external 
radiation (figure 2 (§3.1.2)). This is because the external radiation dissociates a little H 2 in the 
early stage of collapse. However, the difference is negligible. 
3.2-4- Parameter dependence of temperature evolution 

To investigate how parameters affect the evolution of temperature, we systematically 
calculate with changing one of three parameters in the parameter space (no, /, ^21)- For 
unchanged parameters, / = 1.25, n = 10cm~ 3 , and J 2 \ = 10 are used. 

Figure 7 shows the evolution of temperature in the cases with various no. It is seen that 
in the cases with no = 10 — 10 2 cm -3 temperature increases adiabatically since H 2 is photodisso- 
ciated. However, with no higher than 10 2 cm -3 , the filamentary clouds shield themselves from 
dissociation photon and cool owing to H2 cooling. 

Figure 8 shows the evolution of temperature in the cases with various /. In all the cases, 
most of H2 is photodissociated in the early phase (n < 10 2 cm -3 ) and temperature increases 
adiabatically. Each line in figure 8 overlaps each other during the early stage of collapse. 
This is because the effect of photodissociation is similar for the same density. The fragment 
mass depends on /. The filamentary clouds with / < 2.5 fragment during the early adiabatic 
evolution (n < 10 2 cm -3 ). In the case with large /, the filamentary clouds collapse to high 
density since virial temperature is large. The filamentary clouds with / > 2.5 form sufficient 
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amount of H2 to cool during collapse even with strong dissociation radiation (J21 = 10). Once 
the filamentary clouds cool, they continue to collapse and reach the high density (n > 10 13 cm~ 3 ) 
before fragmentation. 

Figure 9 shows the evolution of temperature in the cases with various J21. In the cases 
with the external radiation (J21 > 1), temperature increases adiabatically since most of H2 is 
photodissociated and H2 cooling is suppressed. On the other hand, in the case without the 
external radiation ( J21 = 0), temperature does not increase adiabatically. 

3.3. Fragment mass 

Figure 10 shows the fragment mass for all the parameters by using contours maps in 
n — f plane. Results for different values of J 2 \ = 0, 1, 6.5, and 10 are presented in different 
diagrams. Solid lines in each diagram of figure 10 represent constant fragment mass. The 
dotted line and the dash-dotted line will be referred in §4. In diagram b) with J 2 i = 1 5 the 
region with large fragment mass (> 1O 4 M ) is found in the range uq < 10 L5 cm~ 3 and / < 1.5. 
This region is clearly as the result of the external radiation since such a region does not exist 
in diagram a) with J21 = 0. With larger J21 in diagrams c) and d), it is seen that the region 
with massive fragment (> 1O 4 M ) becomes larger in no — / plane. In diagram c), this region 
spreads up to no ~ 10 2 cm -3 and / = 2.25. In diagram d), this region spreads up to / = 2.5. 

However, in the case with no > 10 2 cm -3 or / > 2.5, it is seen that the fragment mass is 
hardly changed by the external radiation. This is explained as follows : as for the cases with 
large no, the filamentary clouds shield themselves from the dissociation radiation from the early 
stage of evolution. As for the cases with large /, as shown in §3.2.1, the filamentary clouds 
continue to collapse up to density high enough to shield themselves from the external radiation 
and form H2 even if H2 is photodissociated in the early stage of collapse. The filamentary cloud 
in the uniform model whose collapse is homologous tends to collapse to higher density than the 
realistic model whose collapse is run-away collapse. Thus, it is probably as the result of our 
choice of the uniform model that the filamentary cloud with a little larger / than moderate 
value collapses to high density in spite of H 2 loss in the early stage of collapse. In §5, we 
compare the result with the modified one- zone model including the effect of run-away collapse. 
Further investigation including spatial variation will be presented in the separate paper. 

3-4- Effect of supersonic initial velocity 

We assumed initial infall velocity at the cloud surface to be same as the sound velocity. 
Here, we comment on the effect of faster initial infall velocity which may be possible under the 
effect of dark matter gravity. Much faster initial velocity can help the filamentary clouds to 
shield themselves from the dissociation radiation due to the rapid evolution of density before 
fragmentation. We have checked this possibility with the model as in figure 5 (§3.2.2) where 
most of H 2 is photodissociated during the early stage of collapse. It is found that the filamentary 
cloud can collapse to form H2 and shield themselves from the dissociation radiation if the 
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initial velocity is five times or lager than five times sound speed. An example of numerical 
result is shown in figure 11 where photodissociation is suppressed at the early stage of collapse 
and the filamentary cloud does not fragment until density becomes large (~ 10 8 cm -3 ). Thus, 
the filamentary clouds with highly supersonic initial infall velocity tend to avoid the effect 
of the external radiation. This is consistent qualitatively with Hasegawa et al. (2009b) who 
investigated the formation of globular clusters as a result of fragmentation in such a way. 

3. 5. Effect of photodissociation on the initial fraction of H 2 

Although initial H 2 fraction is assumed to be 10 -4 at n = n > 10cm -3 in §3 (and §5), 
this value of H 2 fraction is expected to be affected by the external radiation before the density 
of the cloud reaches n . In previous papers without the external radiation, / H2 = 10~ 4 — 10 -3 is 
adopted as initial states after virialization. In this sense, the above initial value of fn 2 = 10 -4 
should be regarded as for the case where the external radiation turns on at the moment when 
the filamentary cloud with n = no forms. On the other hand, Omukai & Yoshii (2003) considered 
the cases with sufficiently low initial density Uq = 0.1cm -3 , in which the external radiation had 
turned on before the filamentary cloud forms. 

Is our assumption that initial /h 2 is se t to be 10 -4 valid when the external radiation 
turns on before the filamentary cloud forms? In this subsection, we investigate how much _/h 2 
is at n = 10cm -3 in the case where the external radiation turns on at nuv which is lower than 
10cm -3 . We calculate the evolution of H 2 fraction by using the cloud with sufficiently low 
initial density no = 0.1cm -3 . Initial fraction of H 2 is set to be zero. In figure 12, H 2 fraction at 
n = 10cm -3 is shown as a function of nyy for the cases with different J 2 i- It is seen that H 2 
fraction at n = 10cm -3 is much different between in the cases with different J 2 i and nuv- To 
set /h 2 = 10 -4 at n = no in §3 (and §5) is valid only when the external radiation turns on at 
the moment when the filamentary cloud forms. 

4. Analytic investigation 

In this section, we analytically investigate the property shown in the numerical results 
in §3. To explain the property of the collapsing filamentary cloud, three criteria are considered 
in the view point whether or not the filamentary cloud can cool during collapse. 

4-1. Cooling criterion 1 : Whether cooling is effective or not 

There is a critical value n a of initial density that determines whether or not H 2 cooling 
dominates adiabatic heating at the start of collapse. Before we consider the effect of the 
external radiation, we derive n a without the external radiation. If initial density exceeds n a , 
temperature increases from the early stage of collapse and the filamentary clouds fragment into 
massive fragments (~ 5OM ) as shown in figure 2 (§3.1.2). If H 2 cooling dominates adiabatic 
heating at the start of collapse, the following inequality is satisfied : 
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-P d - l -<^. (21) 

dtp p y ! 

The rate of H2 cooling is approximated as 



2.5 x 1(T 2 S 2 /h 2 (-jj^) n « 10 4 cm- 
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Ah 2 - < >» (22) 

8.0 x 1(T 2 Vh 2 (t^J « > 10 4 cm- 3 

(Galli & Palla 1998), where Ah 2 is in units of erg cm _3 s _1 . Assuming that the timescale of 
collapse is the free-fall time (1/p- dp/dt = — tg 1 ), equation (21) becomes 
kbTq 1 A H2 



2nGp < -^ (23) 

pm H " p 

(see Appendix 3). Using Ah 2 for n <C 10 4 cm -3 in equation (22) with /h 2 = 10~ 4 for the initial 

state, the condition for cooling is found to be 

To Y-V/h^ 



n>n a = 1.9xl0 2 cm- J ^77 H^r • (24) 

V3ook; vio-v v ; 

On the other hand, in the case with n 3> 10 4 cm~ 3 , adiabatic heating always dominates H 2 
cooling. Thus, H2 cooling dominates adiabatic heating at the start of collapse for n a < no < 
10 4 cm -3 . In figure 2, since the filamentary cloud has higher initial density than 10 4 cm -4 , 
temperature increases at the early stage of the collapse. Since this condition does not include 
the effect of the external radiation, equation (24) should be accepted as a necessary condition 
for cooling. 

4-2. Equilibrium fraction oj H 2 

Let us prepare to investigate the condition whether H 2 cooling dominates adiabatic 

heating in the early stage with the effect of the external radiation. Since cooling rate depends 

on H 2 fraction, we first estimate the equilibrium fraction of H 2 which is attained when formation 

and photodissociation of H 2 balance under the external radiation. Assuming the chemical 

equilibrium between formation and photodissociation of H 2 , the fraction of H 2 is found to be 

, nf e k u - . . 

/h 2 = -r , (25) 

^2step 

where &h- = 1.0 x 10 -18 T cm 3 s _1 is the reaction rate for H~ channel, 

H + in^H 2 +e-. (26) 

At the initial state in our model, timescale of formation of H 2 , tf or m, is given by 

*for m = ^-r- = 3.33 x 10 18 s (J—) _1 (-^) ' YrferV 1 . (27) 



k H -n f e V300K/ VlOcm- 3 / V10" 4 , 

On the other hand, assuming that Ah 2 is larger than 10 14 cm -2 , timescale of photodissociation 
is given by 
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Equilibrium H 2 fraction, fn 2 ,eq, can be estimated by the condition tf orm = tdh 
Substituting the column density 

M * f P ^ V2 / 1/2 /h 2 

n 
10cm~ 3 




2.40 x 10 2 W 2 (-^y7^yV 2 /H 2 (29) 



into equation (28) 7 , we have 



/H 2 , eq = min 



??. 



, 11/2 / T \ n / 2 

2.88 xlO" 5 ^^- ( — ) J^f'\ 

1 10cm -3 / V300K/ 21 ^ ' 



(30) 



According to equation (30), /H 2 ,cq is expected to be large for high n. In the cases with J 2 i = 
1, 6.5, and 10, equation (30) predicts /H 2 ,eq ~ 1 for n > 76.6cm -3 , 299cm -3 , and 409cm -3 , 
respectively. Especially, /H 2 , eq is expected to be large enough in n> 10 4 cm -3 where H 2 is 
hardly dissociated by the external radiation. In order to check the applicability of equation 
(30) to analytic criteria for massive fragment formation, we compared /h 2 by equation (30) 
with the numerical results of the evolution of the filamentary cloud. In the cases with J 2 i = 6.5 
and 10, it is found that fu 2 ,cq agrees with numerical results within error of 40 %. On the other 
hand, in the case with J 2 i = 1 where H 2 formation dominates photodissociation, it is found that 
/h 2 given by equation (30) is about 2.5 orders of magnitude smaller than the numerical result. 

4-3. Cooling criterion 2 : Whether cooling becomes effective when formation and photodisso- 
ciation of H 2 balance 

In this subsection, we derive the condition whether H 2 cooling dominates adiabatic 

heating in the early stage under the external radiation by assuming that the formation of H 2 

balances with photodissociation. Cooling time is estimated as 

_ 3nk B T 

tcooi - ~7r: ■ {oi) 

2A Ha 

On the other hand, free-fall time of the uniform filamentary cloud is given by 

<» = OT (32) 

By equating t coo i and ts for n <^C 10 4 cm -3 with assuming with /h 2 < 1, we have critical initial 
density n& as 



n ^™-\mK) lit J I' 1 "- < 33 > 

In the case with no > rib, H 2 photodissociation is too weak to halt H 2 cooling. On the other 
hand, the case with n < n& has a possibility to halt H 2 cooling. In diagrams b), c), and d) 

7 There are cases where /h 2 ~ 10 -8 and A^ H2 ~ 2.40 x 10 13 cm -2 < 10 14 cm -2 . However, in such cases, although 
we use / s h = (-/VH 2 /10 14 cm -2 )~ 3 / 4 , we do not face significant error. 



15 



of figure 10, prediction by equation (33) is plotted by the dash-dotted lines 8 . Comparing with 
numerical results, it is found that the dash-dotted line in figure 10 gives us a reasonable criterion 
for above two thermal evolutions. 

4-4- Cooling criterion 3 : Whether cooling is effective with temperature adiabatically increasing 

We investigate another criterion by which H 2 cooling becomes effective during the col- 
lapse under the condition where temperature increases adiabatically as a result of strong pho- 
todissociation of H 2 in the early stage. There are numerical examples presented in figure 4 
(§3.2.1)/figure 5 (§3.2.2) where H 2 cooling is effective/ineffective. Since in both examples most 
of H 2 is once photodissociated, the difference between these two examples seems to be origi- 
nated from the difference in line mass. Here, we derive the critical line mass f c for effective 
H 2 cooling after strong photodissociation from the condition by which the filamentary cloud 
continues to collapse up to the density high enough to shield themselves from the dissociation 
photon. 

We assume n — rib^ 10 4 cm -3 (i.e., formation and photodissociation of H 2 balance). We 
define that H 2 cooling is "effective" if H 2 cooling dominates adiabatic heating when the right 
hand side of equation (5) equals zero (i.e., gravitational force balances with pressure gradient 
force). In the adiabatic evolution, temperature at the density n is represented as 

/ n \ 2/3 

T = T {-) . (34) 

When the right hand side of equation (5) equals zero, we have 

i^ = /(=/ W r )). (35) 

Since the filamentary clouds mainly consist of hydrogen atom, we assume yU~ 1. Using equations 
(34) and (35), we have 

/=£■ (36) 

Using equations (33), (34), and (36), we have 

For the cases with / < f c , H 2 cooling never dominates adiabatic heating and the filamentary 
clouds fragment into very massive fragments (~ 10 4-5 M Q ; c.f., figure 5 (§3.2.2)). In figure 
10, the condition / = f c is shown by dashed lines. It is seen that the dashed line in diagrams 
c) and d) approximately coincides with the solid line for fragmentation mass ~ 10 5 M Q given 
by numerical results of the collapsing filamentary cloud. Thus, we conclude that the condition 
f <f c with equation (37) provides a useful criterion for the formation of very massive fragments. 

8 The dash-dotted line is not drown in diagram a) of figure 10 since we are interested only in the case with 
the external radiation. 
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Assuming equation (30) in the case with J21 = 1, the criterion f c has about factor 3 of error. 

5. Effect of rarefaction wave 

So far, we have assumed a uniform filamentary cloud where density of cloud is constant. 
However, the fragment mass predicted by the uniform model of the filamentary cloud tends 
to be lower than the result with more realistic treatment such as one-dimensional calculation 
(Uehara & Inutsuka 2000). Indeed, the fragment mass in §3 for the cases without the external 
radiation is different from that of Nakamura & Umemura (2002). During collapse, it is expected 
for the density profile to become core-envelope structure with uniform core and rarefied enve- 
lope. One physical reason to form such a density profile is the property of self-gravity Central 
dense region collapses in a shorter time than outer less dense region, and density contrast in- 
creases. This is important for a cold collapsing cloud without pressure. While with the effect of 
pressure, run- away collapse is enhanced even if the initial cloud is uniform. Even without den- 
sity perturbations, pressure gradient force erodes the filamentary cloud from the surface during 
collapse. Rarefaction wave propagates from the outer boundary to the center. Region outside 
the rarefaction wave front is delayed to collapse by the effect of outward pressure gradient force. 
This becomes important in a cloud with non-negligible pressure. In this section, to capture the 
effect of run-away collapse partly, the effect of run-away collapse which is induced by rarefaction 
wave is taken into account in the one-zone model of the collapsing filamentary cloud. Similar 
approach is adopted in a rotating isothermal cloud and is shown to be effective (Tsuribe & 
Inutsuka 1999). To include full characteristics of run-away collapse, the one-dimensional hy- 
drodynamical calculations are required. Results of series of the one-dimensional calculations 
will be reported elsewhere. 

5.1. Modification to the model 

Consider a collapsing filamentary cloud with uniform initial density and pressure. As 
the cloud collapses, radius of the cloud decreases. In addition to this, the rarefaction wave 
propagates inward from the outer boundary according to 

- = -2Trrpc s , (38) 

where f is the position of the rarefaction front, p and / is the line mass and density inside 
f. Combining with the solution of density, velocity for homologous collapse, and I = irr 2 p, we 
can calculate the evolution of f and I. Using f and I, we define and solve the modified virial 
equation instead of equation (5) as 

§ = ~{l-lc{T)}, (39) 

at r 

where v is the infall velocity at the rarefaction wave front. Hereafter we denote this model as 
"rarefied filament model" . Different from the uniform model in previous sections, in the rarefied 
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filament model / decreases as the cloud collapses. Thus, in this model the right hand side of 
equation (39) becomes positive at lower density than the uniform model. Similarly, fragmen- 
tation occurs at lower density in the rarefied filament model. These differences are originated 
from the property of run-away collapse. The same condition as in §2.3 for fragmentation is 
assumed. Mass of the fragment is calculated using Mf rag = 2irrl instead of equation (18). 

5.2. Results 

In figure 13, the result for the rarefied filament model is shown for the case with large line 
mass and low initial density, (f,n ,J 2 i) = (3, 10cm -3 , 0). Compared with the uniform model 
in figure 1, the cloud fragments at lower density as expected. Density at fragmentation is 
5.0 x 10 -12 times that in the uniform model and effective radius of the filamentary cloud at the 
moment of fragmentation is larger by 5.1 x 10 5 . On the other hand, line mass is / = 2.1 x 10~ 2 l 
at the moment of fragmentation. As a result of combination of these effects, the fragmentation 
mass (127M Q ) is 1.1 x 10 3 time larger than in the uniform model. It should be noted that even 
with the rarefied filament model complete property of run-away collapse can not be captured 
for the filamentary cloud with initial density profile with central concentration. In this sense, 
the result of fragment mass still differs by about factor 4 from the previous result of the one- 
dimensional calculation (e.g., Nakamura & Umemura 2002) 9 . 

In figure 14, the result is shown for the case with small line mass and high initial density, 
(/, no, J21) = (1.25, 10 6 cm _3 ,0). Compared with the uniform model in figure 2, in the rarefied 
filament model, fragmentation density is 3.2 x 10~ 2 of that in the uniform model. Effective 
radius of the filamentary cloud at the moment of fragmentation is 5.6 times larger, but line 
mass at that moment is 0.25/. As a result, the fragment mass (64M Q ) is 1.4 times larger than 
in the uniform model. 

Like figure 3, figure 15 shows the fragment mass for various Uq and / in the case with- 
out the external radiation. All lines are similar to each other and can be approximated as 
-Wfrag ~ 30000n ~°' 3 /'~ 4 with an error at most factor 2 at / = 3. Although the fragment mass is 
determined mainly by /, dependence on n$ is stronger than the case in figure 3. This tendency 
agrees with the result of Nakamura & Umemura (2002). 

In figure 16, the result is shown for the case with small line mass, high initial density, and 
strong external radiation, (/, no, J21) = (1.25,10cm -3 , 10). Compared with the uniform model 
in figure 5, in the rarefied filament model, fragmentation density is 0.65 of that in the uniform 
model. Effective radius of the filamentary cloud at the moment of fragmentation is 0.8 of the 
uniform filamentary cloud and the line mass at that moment is 0.29/. As a result, the fragment 
mass is 0.23 of that in the uniform model. Thus, the fragment mass is smaller (~ 2.6 x 1O 4 M ) 
due to smaller line mass, but still larger than without the external radiation. In all of above 
results, thermal evolution is qualitatively similar to that of the uniform model. 



9 Note also that fragmentation timescale of Nakamura & Umemura (2002) is longer by factor 2.5 than ours. 
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5.3. Fragment mass 

In figure 17, contours maps of the fragment mass in no — f plane are shown for the 
cases with J 21 = 0, 1, 6.5, and 10. Solid lines in each diagram represent constant fragment 
mass. Comparing diagram b), c), and d) with diagram a), it is seen that the fragment mass 
for the cases with low initial density (no < 10 2_2 ' 5 cm~ 3 ) is strongly affected by the external 
radiation. In the cases with these low no, fragment mass is mainly determined by no instead of 
/ and massive fragments form in the cases with J 21 > 1. In the cases with high initial density 
(no > 10 3 cm -3 ), the fragment mass is approximately independent of the external radiation. 
Since density is high enough for the filamentary cloud to shield itself from the dissociation 
photon from the early stage of collapse, role of photodissociation is less important. 

We comment on the effects of run-away collapse. By comparing figure 17 with figure 10, 
these effects is clearly noticed. The most remarkable difference is that strong dependence on 
/ in figure 10 becomes weaker in figure 17. Furthermore, sub-solar fragments seen in figure 10 
are not found in figure 17. Thus, we suspect that too small mass of fragments in figure 10 are 
the result of too idealized modeling with the uniform filamentary cloud in previous sections. In 
both of figure 17 and figure 10 massive fragments are seen in the cases with low initial density 
(no < 10 2 cm~ 3 ) and J 2 i (> 1). Formation of massive fragment with the external radiation can 
be regarded as a robust result. 

6. Summary and discussion 

6. 1 . Summary 

In this paper, we investigated collapse and fragmentation of primordial filamentary 
clouds under the external radiation with one-zone models. We numerically calculated the 
thermal and dynamical evolution of the filamentary clouds and estimated the mass of frag- 
ments for a variety of parameters such as rio, f, and J 2 i. According to the uniform model it is 
found that with initial H 2 fraction /h 2 = 10 -4 , low initial density (n < 10 2 cm~ 3 ), and moderate 
line mass (/ < 2) the filamentary cloud loses its cooling ability as a result of photodissociation 
of H 2 by the external radiation whose mean intensity is J 2 i > 6.5. In such a case, gravita- 
tional collapse proceeds adiabatically, and the filamentary clouds fragment into more massive 
fragments (~ 10 4_5 M Q ) than the case without the external radiation (~ 1 — 50M Q ). In the 
cases with lower intensity of the external radiation, the filamentary cloud collapses without 
fragmentation to density which is high enough for H 2 to form as a result of self-shielding. In 
this case, mass of fragments is expected to be similar to the case without the external radiation. 
If the initial density is high (n > 10 2 cm -3 ), the filamentary clouds with moderate line mass 
shields themselves from the dissociation photons. However, in such a high initial density case, 
adiabatic heating dominates cooling. As a result, they fragment into more massive fragments 
(~ 100M Q ) than the low initial density cases with effective H 2 cooling. Summarizing the results 
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of numerical calculations, figure 10 clearly shows that the effect of the external dissociation ra- 
diation increases the fragment mass in low initial density cases (see §3.3). In §4, we derived an 
analytic criterion for the formation of very massive fragments via photodissociation. It is found 
that massive fragment is expected if the cooling time with equilibrium H2 fraction is longer 
than the free-fall time at the end of hypothetical adiabatic collapse. 

In order to modify an unrealistic property of the uniform model where collapse is not 
suppressed until too high density (e.g., ~ 10 15 cm -3 in figure 1), we developed a modified 
version of a simple one-zone model which can partly capture the effect of run-away collapse 
by focusing on the central dense spindle inside the rarefaction wave which comes from the 
outer boundary. According to this rarefied filament model, fragmentation is expected to occur 
before the rarefaction wave front arrives at the center. As a result, fragmentation density is 
smaller and fragment mass is usually larger than in the uniform model. With this new one- 
zone model, we can calculate easily the evolution of a filamentary cloud with run-away collapse 
without assuming free-fall collapse. Different from the uniform model, dependence of fragment 
mass on / becomes weaker, and fragment mass itself becomes larger in the most cases in the 
uq — f plane. This dependence and the value of fragment mass are similar to the result of the 
one-dimensional model (e.g., figure 5a of Nakamura & Umemura 2002). In the uniform model 
and the rarefied filament model, very massive fragments (> 1O 4 M ) form from the filamentary 
clouds with low initial density (n < 10 2 cm -3 ) and the external radiation. This formation of 
very massive fragments is the robust result. 

6.2. Discussion 

We compare the difference of the effect of photodissociation between the spherical cloud 
(Omukai 2001) and the filamentary cloud with radius R. For the uniform spherical cloud, we 
have Ah 2 = n^R oc % 2 n _1 ' 3 oc n 2 ^ 3 and for the uniform filamentary cloud we have Ah 2 = 
(7r/2)nu 2 R oc n^ 2 rT x l 2 oc n 1 / 2 (see §2.2). Since the photodissociation reaction rate festep is 
proportional to N# , we have 

I rT x l 2 sphere 
I n~ 3 / 8 filament. 

Difference in power index in both cases is small and the evolution of temperature would be 
similar to each other. 

We comment on the further evolution of very massive fragments (~ 1O 4_5 M ) which 
form as a result of fragmentation of the filamentary cloud under the external radiation. Since 
each fragment is expected to be nearly spherically symmetric, evolution of spherical cloud under 
the external dissociation radiation will be useful to discuss further evolution of each fragment. 
Susa (2007) investigated collapse of spherical cloud (~ 1O 5 M ) under the UV radiation. When 
the distance between the cloud and a single light source (12OM star) is longer than lOOpc, and 
the light source is turned on when the density of the cloud is 10 2 cm -3 , the author showed that 
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the cloud collapses to 10 6_7 cm -3 . Although our model adopts the uniform external radiation, 
similar evolution may be possible if the radiation field can be regarded as the effect from 
many light sources. In this case, mean intensity depends on the mean distance between the 
filamentary cloud and light sources. In our model with J21 = 6.5, the distance between the 
filamentary cloud in a halo and surrounding star-forming halos (light sources) whose luminosity 
is ~ 10 25 erg s _1 is expected to be longer than 20kpc according to Dijkstra et al. (2008). Thus, 
mean intensity at the surface of the filamentary cloud is expected to be weaker than that of the 
situation of Susa (2007), and the external radiation is not expected to photodissociate H 2 once 
it collapses to high density (~ 10 3_4 cm -3 ). Such a fragment may continue to collapse up to 
10 6_7 cm -3 and shields itself from the external radiation. After that, although H2 will form in 
the collapsing clump and it cools the clump, other effect than cooling physics such as rotation 
and disk formation/accretion may be important. Since the possibility of further fragmentation 
will depends on these processes, the final outcome of each clump is out scope of the present 
paper. Assuming that each clump does not fragment further after it fragment at around the 
loitering point, Omukai & Yoshii (2003) discussed the initial mass function. 

In the realistic situation, the effect of photodissociation is expected to be dominated 
by the single nearest point source (Susa et al. 2009; Hasegawa et al. 2009a). In this paper, 
however, we assume isotropic and steady external radiation field for simplicity. If this source 
of the external radiation is regarded as the group of sources (Dijkstra et al. 2008), new stars 
must be formed continuously around the filamentary cloud. Even in such a case, intensity and 
spectrum may be more complex and evolve with time. Furthermore, the effect from ionization 
photon may not be neglected. These issues may be important but are out scope of this paper. 

Although the rarefied filament model is developed in this paper, it is still questionable 
whether or not all of the filamentary cloud fragment before the moment when rarefaction wave 
reaches the center. To clarify this point, more accurate calculation at least one-dimensional 
hydrodynamical calculation is required. 

We thank Fumio Takahara for fruitful discussion and continuous encouragement, 
Kazuyuki Omukai for showing detailed technical treatment used in Omukai (2001), Shu-ichiro 
Inutsuka for discussion about the rarefied filament model. We also acknowledge the referee for 
improving the manuscript. 

Appendix 1. The opacity of species including deuterium 

According to Omukai (2001), opacity is determined by cross section, partition function, 
and reduced mass. Cross section is the same one that we consider in photoreactions. Reduced 
mass is estimated easily. According to Bron et al. (1973), the ratio of the partition function of 
HD to that of H 2 is estimated. About HD + , the partition function Z is 

Z = Z trans Z TO tZ V ibZ e \ e , (Al) 
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S(x) 



(A2) 



where Z trans oc m 3 / 2 , Z rot oc / oc m, Z vi b ex: m, and Z c i oc m°. Hence Zhd+ = ( m HD+/' m H + ) 7/ ' 2 '^H + - 
Appendix 2. The integral value of equation (14) 

The integrated value of equation (14) can be approximated by 

4vr x<0.01 
-4.10915x 3 + 12.579x 2 - 17.5787a; + 12.5392 0.01 < x < 1 

3.48330X- - 398517 *- - 884854 1< x < 5 

7.57551a;- a264805 *- 1 - 93051 5 < x < 10 

18.3238a;- a236295a; - a258468 10 < x < 20 

48.3927X- - 213846 *- 3 - 28825 20 < x < 30 

54.4202x- a202620a; - 3 - 61369 30 < x < 40 

54.8962x-°' 201651x - a361636 40 < x 

where x = fcjy-R. Error is smaller than 10%. 

Appendix 3. Free-fall time of the uniform filamentary cloud 

Neglecting the effect of kinetic energy and pressure in equation (1) with 

1= f pr 2 dV = -IR 2 , (A3) 

we have 

l^(llR 2 )= -Gl 2 (A4) 

2dt 2 \2 J K J 

^R 2 = -4GI, (A5) 

where / is constant. From equation (A5), R 2 is the quadratic function of t. Hence we can 
express R 2 = at 2 + bt + c (a, b, and c are constant). When t equals zero, R is the initial radius 
and c = -Rq. We assume that the initial velocity is zero and b = 0. According to equation (A5), 



a = —2GI. Using the relation R = Jl/irp, we have 



R=^/-2Glt 2 + l/np. (A6) 

Hence the free-fall time for the filamentary cloud is 

t= vm- (A7) 
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Number 


Name 


Process 


Cross Section (cm ) 


Reference 


al 


H bound-free 


H(n) + 7 -> H+ + e 


7.909 x W- ls n(v/v n )- 3 ;hv n = 13.598eV/n 2 


1 


a2 


He bound-free 


He + 7 -> Hc+ + e 


7.83 x W w [1.66(u/v T )-' 2m - 0.66(i//vr)- 3 - 05 ];hi^ = 24.586eV 


2 


a3 


H~ bound-free 


H" +7^H + e 


10" 18 A ;! (1/A- l/A ) : '/ 2 /(A), A = 1.6419pm, /(A) from cquation(5) of reference 


3 


a4 


Hj bound-free 


Hj+7^H + H+ 


see tablc2 of reference 


4 


a5 


H~ free-free 


H + e + 7^ H + e 


kf(T)kBTn e : kf from equation(6) of reference 


3 


a6 


H free-free 


H++e + 7^H + +e 


3.692 x V,fg s (v,T)v'' i T' 1 l' 2 n c ;'NC take g«(v,T) = 1 




a7 


H 2 -H 2 CIA 


H 2 (i>, J) + H 2 +7 -4 H 2 (o', ./') + H 2 


see figure 1 of reference 


5 


a8 


H 2 -Hc CIA 


H 2 (i>, J) +He + 7-> H 2 (o', ./') + He 


see figurc2 of reference 


5 


a9 


D bound-free 


D + 7^D+ + e 


same as al 


6 


all) 


HD + bound-free 


HD++7^D + H+ 


same as a4 


6 


all 


HD + bound-free 


HD++7^D+ + H 


same as a4 


6 


al2 


D~ bound- free 


D" +7^D + e 


same as a3 


7 


si 


H Raylcigh 


H + 7 -> H + 7' 


5.799 x 10- 29 A- 4 + 1.422 x 10- ;io A- 6 + 2.784 x 10- 32 A" S 


8 


s2 


Tomson 


e + 7 — > e + 7' 


6.65 x 1CT 25 


1 



Table 1. Continuum processes. The wave length A is in units of 10~ 6 cm. REFERENCES-(l)Rybicki & Lightman 1979; (2) 
Osterbrock 1989; (3) John 1988; (4) Stancil 1994; (5) Borysow, Jorgensen, & Zheng 1997; (6) Galli & Palla 1998; (7) Frolov 
2004; (8) Kurucz 1970. al-a8, si, and s2 are considered in Omukai 2001. 
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Fig. 1. Evolution of the temperature (top), /h 2 (the second from top), I and l c (§2.1:third), the 
heating and cooling rate (fourth), and tdyn and tf rag (§2.3:bottom), respectively, as a function of 
the density for model with (/, n , J21) = (3, 10cm -3 , 0), in which "adiabatic" denotes the adia- 
batic heating, "H2" does the H2 line cooling, "HD" docs the HD line cooling, and "chemical" 
does the chemical heating or cooling. We omit the continuum cooling because it is not effective. 
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Fig. 2. Same as figurel, but (/,n , J 21) = (1.25,10 6 cm- 3 ,0). 



&. 10' 




Fig. 3. The fragment mass in the cases without the external radiation with various of no and /. The 
hollowed region in the cases with / = 2.5 and no > 10 5 cm~ 3 in which H2 forms by three body reaction cools 
strongly. Since cooling helps collapse, the filamentary clouds collapse up to high density (~ 10 13 cm~ 3 ). 
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Fig. 4. Same as figurcl, but (f,no, J21) = (3,10cm 3 ,10). 
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Fig. 5. Same as figurcl, but (f,no, J21) = (1.25,10cm 3 ,10). 
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Fig. 6. Same as figurel, but (f,n , J 21) = (1.25,10 6 cm- 3 ,10). 
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Fig. 7. The evolution of the temperature for the cases with the variety of initial den- 
sity, fixed /, and fixed J21. The parentheses in the figure denotes the parame- 
ters (/, ng, J21). The right break in each line denotes (T, n) at the fragmentation. 
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Fig. 8. Same as figure7, but for the cases with the variety of line mass parameter, 
fixed no, and fixed 3i\. Cross denotes (T, n) at the fragmentation. It is clearly seen 
that fragmentation is divided into two groups with n ~ 10 2 cm~ 3 and n > 10 13 cm -3 . 
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Fig. 9. Same as figurc7, but for the cases with the variety of intensity of the external ra- 
diation, fixed no, and fixed /. The lines for J21 = 6.5 and 10 overlap each other. 
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Fig. 10. The contours map for the fragment mass in no — / plane for the case with a) J%\ = 0, 
b) J21 = 1, c) J21 = 6.5, and d) J21 = 10. The number near each solid line is mass of 
fragment in units of Mq. The dashed line and the dot-dashed line denote equation (38) 
and equation (34), respectively. In the region on right of the dash-dotted line, i coo i > iff. 



32 



10° 



10 2 
10 -3 

J 10* 

101 

10^ 

10 a 
1CT 8 



,10' 



10^ 
10° 



10" - 



10" 4 1 



10;: 
1fl 



'I I I 'I 




H'V 




-J I l_d 1 ■-J I . .1 



n(cm -3 ) 



Fig. 11. Same as figurel, but the initial velocity is five times of the sound speed and 
(/, n , J 21 ) (1.25, 10cm- 3 , 10). 
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Fig. 12. The fraction of H2 at n = 10cm 3 , as function of nuv which is the 
density at which light sources turn on. Each line corresponds to various J 2 i- 
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Fig. 13. Same as figurel, but for the rarefied filament model, (/,no,«/2i) = (3,10cm 3 ,0). 
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Fig. 14. Same as figure2, but for the rarefied filament model, (/,no, J21) = (1.25,10 6 cm 3 ,0). 
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Fig. 15. The fragment mass in the cases without the external ra- 

diation with various of rig and / in rarefied filament model. 



35 



itr 



10 -3 

10"! 

10 ^ 

10~ 7 



10" 




10 



v 10' 

to 

CO 

I 10* 

01 ,1 

.a. 

< 10" 5 



10 b 

10 10 

c 10 9 

§ 10 8 

10 7 



10° 



10 1 



■c 

I 

adiabatic 
H 2 
HD 

chemical 

Urn 
T frag 



n(cm 3 ) 



10' 



Fig. 16. Same as figure5, but for the rarefied filament model, (f,no, J21) = (1.25, f 0cm 3 ,10). 
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Fig. 17. The contours map for the fragment mass for the case with a) J21 — 0, b) J21 = 1, c) 
J21 = 6.5, and d) J21 = 10. The number near each solid line is mass of fragment in units of Mq. 
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